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Optimal working conditions for thermoelectric generators with 
realistic thermal coupling 
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PACS 84 . 60 . Rb - Thermoelectric energy conversion 
PACS 85 . 80 . Fi - Thermoelectric devices 
PACS 88 . 05 . De - Thermodynamic constraints 

Abstract -We study how maximum output power can be obtained from a thermoelectric gener- 
ator (TEG) with nonideal heat exchangers. We demonstrate with an analytic approach based on 
a force-flux formalism that the sole improvement of the intrinsic characteristics of thermoelectric 
modules including the enhancement of the figure of merit is of limited interest: the constraints 
imposed by the working conditions of the TEG must be considered on the same footing. Introduc- 
ing an effective thermal conductance we derive the conditions which permit maximization of both 
efficiency and power production of the TEG dissipatively coupled to heat reservoirs. Thermal 
impedance matching must be accounted for as well as electrical impedance matching in order to 
maximize the output power. Our calculations also show that the thermal impedance does not only 
depend on the thermal conductivity at zero electrical current: it also depends on the TEG figure 
of merit. Our analysis thus yields both electrical and thermal conditions permitting optimal use 
of a thermoelectric generator. 



Introduction. — Thermoelectric power generation 
poses challenges which are of fundamental and technolog- 
ical nature. The development of efficient thermoelectric 
systems is widely recognized as a strategic topic of applied 
research in view of problems related to, e.g., waste heat re- 
covery and conversion to electricity. This is reflected by an 
abundant litterature on the subject pQ. Recent progress 
in technological development of TEGs has relied on ad- 
vances in material sciences: new materials and new tech- 
niques to produce specific structures have permitted the 
improvement of device performance through the charac- 
terization and optimization of their electrical and thermal 
transport properties (see the review of Di Salvio [2] and 
the recent one of Shakouri [3] ) . The performances still are 
quite modest though and, as for all heat engines, thermo- 
electric generators are subjected to the laws of thermody- 
namics, which impose an upper bound to their efficiency, 
the so-called Carnot efficiency, rjc- Much effort thus is 
invested to seek ways to improve the intrinsic properties 
of thermoelectric modules and hence approach the upper 



efficiency limit. 

These properties often are summarized into the paradig- 
matic figure of merit ZT ', which characterizes the perfor- 
mance of a device at average temperature T [1] (a pre- 
cise definition of the quantity Z is given in the next sec- 
tion). As an illustration of progress we mention the recent 
results of Snyder and coworkers: with p-type doping of 
PbTe-based semiconductor materials, i.e. through band 
structure engineering, a value of ZT as high as 1.8 has 
been reported [4]. This impressive value was obtained at 
a temperature of 850 K, which means that much progress 
remains to be done to achieve this level of performance 
at room temperature. Another aspect of the problem 
was pointed out by Nemir and Beck [5]: since the fig- 
ure of merit is determined by a set of three parameters 
(the electrical conductivity, the thermal conductivity, and 
the thermopower), there is in principle an infinite num- 
ber of possibilities to obtain a given value of ZT, and one 
direct consequence of this fact is that ZT alone is insuf- 
ficient to characterize the performance of a thermoelec- 
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trie system. Furthermore, as discussed below, the abil- 
ity for a TEG to deliver a high output power also de- 
pends on impedance matching which must account for 
thermal resistance; this point directly relates to the more 
general problem of the optimization of the working condi- 
tions of a non-endoreversible thermodynamic engine under 
specific constraints, which indeed differ from those of en- 
dorcvcrsiblc engines, also known as the Novikov-Curzon- 
Ahlborn (NCA) conditions [BHS]. 

In this work, we thus adopt a particular viewpoint: we 
do not seek ways to obtain ever higher values of ZT; 
rather, we want to understand how one can optimize both 
electrical and thermal conditions so that a thcrmogener- 
ator may produce a maximum ouput power. Indeed, if 
the connections of a TEG to two heat reservoirs are as- 
sumed to be ideal, the best working conditions are well 
known; but when thermal dissipative couplings exist, as 
for all real systems, these conditions turn out to be dif- 
ferent. A global reflexion on realistic working conditions 
of TEGs thus is necessary to advantageously exploit the 
benefits of the works on the materials properties. Much 
work has already been devoted to the modeling and sim- 
ulation of thermoelectric devices, but the thermodynamic 
optimization has been a subject of debate, especially when 
considering the maximal efficiency versus maximal power 
strategy [TUHE]. 

In our work, we consider that the electrons in the TEG 
form a carrier gas that can be described in the frame of 
linear out-of-equilibrium thermodynamics developped by 
Onsagcr [13 /13] and Callen [15]ll6j . This approach is most 
convenient to study on the same footing the thermody- 
namic forces inducing the irreversible processes in a sys- 
tem, and the response of this system in terms of fluxes. 
If the fluctuations are sufficiently small a linear force-flux 
coupling provides a correct description of a TEG since the 
physics of such systems is based on the interplay between 
Ohm's law and Fourier's law. Our article is organized 
as follows. In section 2, we introduce the model ther- 
moelectric generator, the definitions and notations we use 
throughout the paper. In Section 3, we analyze the electri- 
cal and thermal conditions for output power and efficiency 
maximization considering non-ideal thermal contacts and 
a fixed figure of merit. We end the paper with a discussion 
and concluding remarks. 

Model of thermoelectric generator. We con- 
sider a single-leg module placed between two ideal heat 
reservoirs as depicted in Fig. Q] (the doping type does not 
affect the generality of the conclusions of our study). The 
temperatures of the reservoirs are denoted T co id (colder 
one) and Thot (hotter one) respectively. The two thermal 
contacts are characterized by the thermal conductances 
-Kcoid and ifhot, respectively. The TEG is characterized by 
its isothermal electrical resistance R, its thermal conduc- 
tance if teg, and the Seebeck coefficient a. The thermal 
conductance -FTteg reduces to the conductance K v=0 , un- 
der zero voltage (electrical short circuit), and to the con- 




Fig. 1: Thermoelectrical (left) and thermodynamical (right) 
pictures of the thermoelectric module and the load. 

ductancc K I=0 , at zero electrical current (open circuit). 
Electrons and phonons contribute to the thermal conduc- 
tance, and Z is given by : Z = a 2 / RK I=Q . 

Thermal and electrical currents. In real systems, the 
heat exchangers between the temperature reservoirs and 
the thermoelectric engine contribute to the energy loss and 
the related decrease of efficiency. In other words, entropy 
is produced in the engine as well as in the exchangers. 
The thermodynamic system we consider is more complex 
than an NCA one, and assuming that the incoming and 
outgoing heat fluxes are linear in the temperature differ- 
ence (see Ref. [T7] for further detail on the laws governing 
heat fluxes), we describe the TEG characteristics with the 
force-flux formalism, which yields the following equation: 

/ I \ 1 / 1 _ a \ / AV \ 

\ Iq J R \ aT a 2 T + RK 1=a ) \ AT' J ' 

(1) 

where / is the electrical current through the load, and Iq 
is the thermal current; AV and AT' = X^m — T c m are the 
voltage and the temperature difference across the TEG, 
ThM and T c m being respectively the temperatures on the 
hot and cold sides of the TEG. The open-circuit voltage 
is V oc = a AT'. The average temperature in the module 
is taken as T = (T c m + T\im)/2. Since we assume that 
the system's response is linear, the temperature difference 
AT' is necessarily small compared to the mean temper- 
ature. We show below how an expression for the TEG 
thermal conductance .Kteg can be derived using two dif- 
ferent ways. 

Within the force-flux formalism the thermal current is 
expressed as the sum of the contributions of convective 
heat transfer and steady-state conduction: 

I Q = aTI + K^ B AT' (2) 

Assuming that the load is simply resistive, Ohm's law ap- 
plies as follows: AV" = —R\ondI and the electrical current 
I reads: 

i _ AV + oAT' _ oAT' ^ 

R -Rload + R 
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This expression of the electrical current / is fed to Eq. ^ , 
which defines the TEG thermal conductance Ateg : 



la = 



a 2 T 



R 



load 



R 



K, 



AT' = AtegAT' (4) 



We see that the thermal conductivity Ateg depends on 
the electrical operating point since its expression given 
above contains the load electrical resistance. 

Now, it is instructive to derive Ateg in a different fash- 
ion. The relationship between the two thermal conduc- 



tances K v=0 and K I=a of the TEG 



K V=0 =K 1=0 (1 + ZT) 



(5) 



can be extended to the following phenomenological for- 
mula: 

K TEG (I) = K I=0 (i + -Lzrj , (6) 

where J gc = aAT'/R, is the short circuit current such 
that: Ateg (Ac) = K v=0 . It is easy to check that equa- 
tions © and Q yield exactly the same expression for 
Ateg- This approach is interesting for two reasons: first, 
the phenomenological law © is validated, not solely be- 
cause the relationship between the two thermal conduc- 
tances K v=a and K I=0 is recovered for I — I sc ; second, it 
is important to note that since the short circuit current I sc 
depends on the effective temperature difference across the 
TEG, AT', there is no closed form solution for the global 
currents and potentials distributions in the TEG. 

Temperatures across the TEG. The analysis devel- 
opped so far assumes the knowledge of AT', but it is 
more useful to obtain expressions of power and efficiency 
as functions of the temperature difference between the two 
reservoirs, AT = Thot — T co id. Here, we give a brief out- 
line of the calculations that yield the relationships between 
T c m, ThM, T co id, and Thot- First, we define the incoming 
heat flux: Q m = Ah t(Th t — ThM) and outgoing heat flux: 
Qout = A co i d (T cM - T cold ); following Ioffe's approach [20], 
these may also be written as: 

Qin = aT hM I - ^RI 2 + A I=0 (T hM - T cM ) (7) 
Qout = aT cM I+^RI 2 + K I=a {T hM -T cM ) (8) 

These equations yield a 2 x 2 system which links T c m and 
ThM to T co id and Thot.: 



V 



hot 



1 RI \ 
2" Ahot 



—T co id — 



I RI 



2 A r 



.1.1 



/ 



M 2 i M22 




(9) 

where the four dimensionless matrix elements are given by: 
M n = A J=n /A ho t + aI/K hot + I, M12 = -K I=0 /K hot , 
M 2 i = A J=0 /A cold , M22 = aI/K cold - A I=0 /A co id - 1- 



The analytic expressions of the temperatures ThM and T c m 
are easily obtained by matrix inversion, but the exact ex- 
pression of AT' as function of Thot and T co id is cumber- 
some, and not necessary for the discussion in the subse- 
quent part of the article. Instead, it is worthwhile to seek 
an approximate but straightforward relationship between 
AT' and AT, which we do as follows. Introducing the total 
contact conductance, A CO ntact, as A~ ntact = K~ ] d + K^ v 
and assuming that the thermal flux is constant in the 
whole system, we obtain the following simple relation be- 
tween AT' and AT: 



AT' = T, 



hM 



T 



A r 



M 



ntact AJ , ^ 

AtEG + Acontact 

using an analogue of the voltage divider formula. The 
assumption we just made amounts to consider that the 
produced electrical power is negligible in comparison to 
the thermal current entering the thermoelectric genera- 
tor, and it thus holds well when the temperature differ- 
ence AT is not too large: in this case the Carnot effi- 
ciency is small, and the real efficiency even much smaller. 
Note that equation (fT0| still holds for the dissymetric case 
(A co id 7^ A'hot). In such situation, the dissymetry is con- 
tained in T: the mean temperature gets closer to Thot 
(Tcoid) if A co id is smaller (greater) than Ahot- However, 
these variations are small since T is comprized between 
T co id and Thot = T co id + AT: T is approximately equal to 
the mean temperature between heat reservoirs. 

Thevenin generator model. In Fig.[T]the electrical part 
of the TEG is viewed as the association of a perfect gener- 
ator and a resistance which is the physical resistance of the 
generator. The open circuit voltage depends on the tem- 
perature difference seen by the TEG: V oc = a AT'. In the 
presence of finite thermal contacts the temperature differ- 
ence AT' depends on the electrical load, hence the tension 
generator can no longer be considered as perfect since its 
characteristics depend on the load. To express this depen- 
dence explicitly, we insert the expression of Ateg given 
by Eq. (J6j> and the definition of the short circuit current 
J sc into Eq. (|I0[) , and we find V oc as the sum of two terms: 



V nc = oAT 



K< 



-IR- 



ZT 



Kl=o + -^contact 1 + A CO ntact/ A /=Q 



(11) 



The first term on the right hand side is independent of 
the electrical load, the second depends on the electrical 
current delivered: V oc = V^' c — IR' ■ We thus obtain a 
rigorous Thevenin modeling of the electrical part of the 
TEG with the definitions of the open circuit voltage given 
by Vq C and the internal resistance is Ateg = R + R'- 

Maximization of power and efficiency with non- 
ideal thermal contacts and fixed ZT. — 

Maximization of power by electrical impedance match- 
ing. The electrical power produced by the TEG can be 
simply expressed as 



P = VI 



V!/R 



load 



(A 



TEG 



R\ 



oad 



|2 ' 



(12) 
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The maximization of the produced output power for a 
given thermal configuration therefore corresponds to: 



1? 



load 



R 



TEG, 



(13) 



which expressed in a more conventional way using the ratio 
m = Rioad/R defined by Ioffe [20], reads: 



ZT 



K con t ac t/K I=Q + 1 



(14) 



We see, as did Freunek and co-workers [21], that the 
electrical impedance matching (|14[) does not correspond 
to the condition m = 1 (or, equivalcntly, i?i oa d = R) of 
the ideal case since the equivalent resistance -Rteg of the 
generator has an additional part due to the finite thermal 
contact coupling. When the electrical resistance matching 
is satisfied, the maximum ouput power reads: 



(-^contact AT) Z 



ZT 



i(K I=0 + K contact )T 1 + ZT + K coatact /K I=0 

(15) 

Maximization of power by thermal impedance matching. 

If we suppose that the TEG is used in a particular en- 
vironment which imposes fixed thermal conductances for 
the contacts, we have to answer the question: how the 
thermal properties of the TEG can be chosen so that a 
maximum output power is obtained? This question di- 
rectly relates to the general problem of the optimization of 
the working conditions of a non-endoreversible engine with 
heat exchangers coupled to the temperature reservoirs. 
This framework extends the classical so-called Novikov- 
Curzon- Ahlborn configuration [6 9 specialized to endore- 
versible engines, which implies that the heat exchangers 
are the only location for entropy production, a process 
thus governed by only one degree of freedom. In a non- 
cndorcvcrsiblc engine, entropy is produced inside the en- 
gine, and so in two different ways: the Joule effect and 
the thermal conduction effect; this confers two additional 
degrees of freedom to the system. 

If we consider that ^contact is fixed by an external con- 
straint, optimization of power may be achieved with re- 
spect to K I=0 ; calculations yield the condition: 



Kc 



= 1 



ZT 



which corresponds to the equality: 



K 



TEG 



(16) 



(17) 



Equation (|17p is similar to that derived by Stevens in 
Rcf. [23] where the thermal impedance matching corre- 
sponds to the equality between the thermal contact resis- 
tance and that of the TEG; however the difference with 
our result above is that the thermal resistance used in 
Ref. [23] for the thermoelectric module is obtained under 
open circuit condition and, as such, it does not account for 
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Fig. 2: Maximum power as function of the ratio K I=0 /K contact 
for various ZT values at fixed ifcontact. In the inset, the curves 
(with ideal and finite thermal contacts) are computed with the 
data of Ref. [5] where the authors studied P max for three values 



of K, 
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the convective part of the thermal current, while -KTteg 
defined in Eq. (0} does. We end this part by highlighting 
the symmetry between electrical and thermal impedance 
matching respectively given by equations (|14p and (|16[) . 

Simultaneous thermal and electrical impedance match- 
ing. The optimal point for a TEG offering the possibil- 
ity to use two degrees of freedom: one electrical, i?i oa d, the 
other thermal, K I=0 (still with ZT fixed), in a particular 
configuration imposed by the environnement, is found by 
joint optimization of the electrical and thermal conditions. 
In other words, equations ([Ft]) and (fTB"]) have to be solved 
simultaneously. We find: 



Kc 



K, 



\J~ZT 



+ 1 



\JZT- 



1 



(18) 
(19) 



We note that Eq. ([T5]) was put forward by Freunek and 
co-workers |21j . and that Yazawa and Shakouri obtained 
both equations [35]. With these two impedance matching 
conditions, wc find that the maximum power produced by 
the TEG is given by: 



P, 



t ZT (AT) 2 



i + vT+zt 



m\ 4T 



(20) 



On the importance of thermal impedance matching. 
The variations of the maximum power P max as a func- 
tion of the ratio K I=g /K contac t [Eq. (|15p] arc shown in 
Fig. [2] for three values of the figure of merit ZT. For a 
given ZT, P max displays a bell-shape appearance (notice 
the use of a logarithmic scale for the abscissa axis). As can 
be expected, higher values of ZT yield greater values for 
the maximum of -Pmax and larger widths at half maximum. 
The maxima are shifted towards the region of lower values 
of the ratio K I=0 /K. contact- Figure [2] demonstrates the im- 
portance of thermal impedance matching: a high value of 
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ZT does not guarantee a greater P max for any value of the 
thermal conductance at zero electrical current K I=0 : for 
instance, P max at K t=0 = K contact for ZT = 1 is greater 



than P max at K, 



5K C 



for ZT = 10. 



In the inset of Fig. [21 two curves represent the maxi- 
mum power as a function of K I=0 for finite and perfect 
thermal contacts respectively; these shed light on the ob- 
servations of Nemir and Beck [5] , which we mentionned in 
the Introduction: to analyze the impact of thermal con- 
tacts on device performance, they considered various con- 
figurations giving the same value for the figure of merit 
ZT. They concluded that for a given value of contact 
thermal conductance the impact on the performance is 
strongly influenced by how the fixed figure of merit of the 
thermoelectric module is achieved. With our analysis one 
can now understand why the TEG with the highest K I=0 
presents the largest performance degradation. 

For fixed values of the electrical resistance and Seebeck 
coefficient of the module, electrical impedance matching 
implies that ZT is not constant, hence the optimization 
amounts to obtain the lowest possible value of the thermal 
conductivity inducing an increase of ZT ', which overcom- 
pensates the mismatching. 

Maximum efficiency. The conversion of the heat cur- 
rent Iq into the electric power P is a process characterized 
by the efficiency rj = P/Iq, which we write as: 



V 



-^contact + KteG P 
i^contact-^TBG AT 



(21) 



considering Eqs. @ and (|10l) . If the thermal contacts are 
ideal, then the expression above reduces to 



1 + m + (ZT hot ) 1 (1 + m) 2 - ?7c/2 



(22) 



We find that the value of m (or equivalently the load 
resistance) which maximizes the efficiency (|21[) is: 



m 




1 + ZT) [l + ZT 



K, 



(23) 



It explicitly depends on the thermal conductances Acontact 
and K I=0 . If the working conditions lead to modifica- 
tions of A'contact (as, e.g., for liquid-gas heat exchangers), 
the operating point of the thermoelectric device changes 
consequently. It is thus interesting to see if m v=llm ^ is 
bounded when the ratio K I=0 / K contact varies. We checked 
from Eq. Q that the mean temperature T varies very little 
with -Kcontact so we may safely consider that the figure of 
merit is fixed without loss of generality for the discussion 
that follows. 

Analysis of optimization and power-efficiency trade-off. 
Accounting for finite thermal contacts in the TEG model 
induces changes in optimal values of the electrical load to 
achieve maximum power or efficiency. The optimal param- 
eters TO pt [ m ,,=, max f° r maximum efficiency in Eq. ()23l) . 
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Fig. 3: Variations of the optimal parameters m (dashed 
line) and m p=p (dashed-dotted line) as functions of -ft'eontact 
scaled to K I=0 , for ZT = 1 and ZT = 3. The shaded areas 
corresponds to the region of best optimization. 



and m P=P for maximum power, in Eq. (|14[) ] arc plotted 
against A con tact (scaled to K I=0 ) in Fig. [3] considering two 
values of the figure of merit: ZT = 1 and ZT = 3. For a 
given ZT and A con tact S> A' J=0 (conditions close to per- 
fect thermal contacts) the maximum power and maximum 
efficiency are well separated: m — > yl-j- ZT, and 
TO p=p max — > 1 (the separation between both obviously 
increases with ZT). 

Conversely, for A contac t < K I=0 we obtain TO„ = „ max — > 
1 + ZT, which also is the upper bound to m p _ p : both 
optimal parameters coincide. At first glance, the conver- 
gence of both optimal parameters towards the same value 
can be seen as valuable since this implies that there is no 
compromise to make between efficiency and power; how- 
ever this is also results in a large performance decrease. 
Actually the regions lying between each pair of curves can 
be considered as the optimal regions to satisfy the power- 
efficiency trade-off. From this point of view we see that the 
narrowing of this zone, which also comes along with lower 
values of ZT, is not at all desirable as it offers less flexi- 
bility in terms of working conditions of the thermoelectric 
generator. 

The figure U] displays two power-efficiency curves, one 



for K c 



^contact = 10A /=Q , tne ouier lor -t\ C ontact — ^^ 7=0 - 

The narrowing of the optimal zone as the contact thermal 
conductance decreases is evidenced and thus confirms the 
observations made for Fig. [31 The arrows indicate the 
maximal values of P max and ^ max in the two cases. The 
maximum power is more sensitive than the maximum ef- 
ficiency to the quality of the thermal contact: in the ex- 
ample we give, the ratio of the two highest values of -Pmax 
(one for each curve) is around 6, while the ratio of the two 
highest values of T) max is around 1.15. When the contact 
is not good enough the power-efficiency curve reduces to a 
point located at the origin where the efficiency and power 
optimisation are indeed identical. 



the other for K c , 
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Efficiency 



Fig. 4: Power versus effcicency curves for two cases with a fixed 
figure of merit ZT = f . Only the value of ^contact varies. 

Discussion and conclusion. Increasing the effi- 
ciency of thermoelectric conversion processes in real de- 
vices remains a topical research problem after several 
decades of efforts and progress. Minimization of entropy 
production methods have proved fruitful along the years 
in various sectors of thermal engineering and sciences [5] 
and new contributions are added to this field with, e.g., 
the recent introduction of the thermoelectric potential 
111. 

Thermoelectric efficiency has been theoretically shown 
to approach the Carnot efficiency in the case of an ergodic 
Lorentz gas with a large figure of merit [25j , but for real en- 
gines it is the efficiency at maximum power that is sought; 
such efficiency was initially shown to be bounded by the 
so-called formula of Curzon and Ahlborn in the specific 
case of an endoreversible heat engine [5] in the frame of 
finite-time themodynamics. The analysis of Curzon and 
Ahlborn was later put on firmer grounds in the frame of 
linear irreversible thermodynamics assuming a strong cou- 
pling between the heat flux and the work ; applications 
to a nanoscopic quantum dot system |27] and extension to 
stochastic heat engines [25] followed. All these works as- 
sume strong or perfect couplings. 

Using a force-flux formalism, we obtained the thermal 
and electrical conditions which allow the maximum power 
production by a thermoelectric generator non-ideally cou- 
pled to heat reservoirs. Introducing an equivalent thermal 
conductance for the generator, which includes conductive 
and convective heat transports, we showed that the ther- 
mal impedance matching can be expressed in a simple 
fashion: the equality of the contact thermal conductances 
and the equivalent thermal conductance. Our analysis is 
thus physically more transparent than that based on the 
geometry of the elements that compose the thermoelectric 
generator. 

Our calculations show that the interplay between the 
thermal and electrical properties of the TEG makes diffi- 
cult the search for the optimum conditions for maximum 
output power production. The importance of the quality 
of the contacts between the TEG and the heat reservoirs 



is demonstrated: high values of ZT are of limited inter- 
est otherwise; in fact, one should search for the ZT that 
allows both electrical and thermal impedance matching, 
which is also the basic idea of the compatibility approach 
|18j for ideal thermoelectric systems. 

* * * 

This work is part of the CERES2 and ISIS projects 
funded by the Agence Nationale de la Recherche. We 
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